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Abstract. In order to characterize flux flow through disordered type-II superconductors, 
we investigate the effects of columnar and point defects on the vortex velocity / voltage 
power spectrum in the driven non-equilibrium steady state. We employ three-dimensional 
Metropolis Monte Carlo simulations to measure relevant physical observables including the 
force-velocity / current-voltage (I-V) characteristics, vortex spatial arrangement and struc- 
ture factor, and mean flux line radius of gyration. Our simulation results compare well to 
earlier flndings and physical intuition. We focus speciflcally on the voltage noise power spec- 
tra in conjunction with the vortex structure factor in the presence of weak columnar and 
point pinning centers. We investigate the vortex washboard noise peak and associated higher 
harmonics, and show that the intensity ratios of the washboard harmonics are determined 
by the strength of the material defects rather than the type of pins present. Through varying 
columnar defect lengths and pinning strengths as well as magnetic flux density we further 
explore the effect of the material defects on vortex transport. It is demonstrated that the 
radius of gyration displays quantitatively unique features that depend characteristically on 
the type of material defects present in the sample. 

PACS. 74.25.Qt Vortex lattices, flux pinning, flux creep - 74.25.Sv Critical currents - 
74.40.-|-k Fluctuations (noise, chaos, nonequilibrium superconductivity, localization, etc.) 

1 Introduction a significant role in their motion above that thresh- 
old, providing, on a coarse-grained description level, 

The physics of interacting vortex lines in high-tem- an effective friction or viscosity, 

perature superconductors subject to strong thermal Aside from the obvious relevance for techno- 

fluctuations and pointlike or extended disorder is logical applications, driven magnetic flux lines in 

amazingly rich and has been a major research fo- type-II superconductors also represent one of the 

cus in condensed matter physics in the past two few cleanly experimentally realizable systems of in- 

decades [T]. Vortex motion plays a crucial role in teracting particles in a non-trivial non-equilibrium 

the transport properties of type-II superconductors steady state. A thorough understanding of the ensu- 

in external magnetic fields. In the presence of a suf- ing phase diagram and full characterization of each 

ficiently large applied current magnetic flux lines emerging steady state should shed light on the rich 

will experience a Lorentz force and drift perpen- and still rather incompletely understood features of 

dicular to both the current and applied magnetic non-equilibrium systems in general. Since stochas- 

field. The motion of these current-encircled mag- tic fluctuations and intrinsic correlations typically 

netic flux filaments induces an electric field parallel play a significant role away from thermal equilib- 

to the applied current resulting in power dissipa- rium, it is desirable for each emerging stationary 

tion and hence an Ohmic voltage drop across the state to attain a thorough quantitative understand- 

material proportional to the velocity of the vortices, ing of the fluctuations. 

While an externally applied current serves as a driv- In superconductors specifically, pinning effects 

ing force, underlying defects in the superconducting on vortex motion are also reflected in the voltage 

material inhibit vortex motion. Material defects pin noise power spectrum [2*. For instance, slightly abo- 

vortices below a critical applied force, and also play ve the critical force, in the presence of strong disor- 
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der, depinning of vortices is observed to proceed via 
flow through plastic channels These 'rivers' of 
vortices form at different locations in the sample, 
and flow around 'islands' of temporarily trapped 
flux lines resulting in incoherent motion. Such be- 
havior has been observed experimentally [3] , as well 
as in two-dimensional computer simulations [5] , and 
is characterized in the velocity or voltage frequency 
power spectrum by a broadband noise signal which 
obeys a power law, as also demonstrated in 

recent three-dimensional numerical work [6]. 

Well above the critical force it has been observed 
that vortices are more translationally ordered than 
at low velocities [7]. It might be expected that at 
a sufficiently high drive, the vortices would form a 
moving Abrikosov lattice since the effective pinning 
force from the disorder on each vortex varies rapidly 
and would therefore be less effective [8j. However, 
it has been shown by Le Doussal and Giamarchi [9] 
that some modes of the disorder are not affected 
by the motion even at large velocities. As a result 
the vortices enter what has been termed a 'mov- 
ing glass' phase. Here, subtle competitions between 
elastic energy, disorder, and dissipation lead to the 
transverse displacements becoming pinned into pre- 
ferred time-independent configurations resulting in 
stationary two-dimensional channels. In the mov- 
ing glass phase, vortices thus follow each other in a 
manner similar to beads on a wire. 

There exist a few possible coupling regimes be- 
tween these elastic channels. For strong point disor- 
der approximately parallel elastic channels are com- 
pletely decoupled, while the periodicity along the 
direction transverse to the drive is maintained. Here 
sharp delta-function Bragg peaks with nonzero reci- 
procal-lattice vector components along the direc- 
tion of motion are lost, while peaks with only trans- 
verse components remain. This regime is known as 
the 'moving transverse glass' or 'moving smectic' 
phase [5] , and is supported by recent numerical sim- 
ulations in two and three dimensions |10|llj . On 
the other hand, for weak point disorder, or large 
velocities, relative deformations grow only logarith- 
mically with distance; hence, the vortex structure 
maintains quasi long-range order corresponding to 
complete elastic coupling between vortices. This state 
is known as a 'moving Bragg glass' and is charac- 
terized by algebraically divergent structure factor 
peaks at small reciprocal lattice vectors [9]. 

The interaction of the moving Bragg glass with 
the underlying material defects is manifested by a 
characteristic peak in the power spectrum corre- 
sponding to the periodicity of the vortex lattice [9] . 
Defects in the material temporarily slow vortices 
resulting in 'stick-slip' motion. In a structurally or- 
dered phase such as the moving Bragg glass, this be- 
havior is repeated, resulting in a periodically vary- 
ing average overall velocity. With the lattice vector 
of the Bragg glass oriented in the flow direction, the 



resulting characteristic frequency associated with 
this motion is known as the 'washboard' frequency. 
Related phenomena are certainly not unique to dri- 
ven flux lines in type-II superconductors, but are, 
for example, well-established in charge- and spin- 
density wave systems [12]. 

Random point disorder need not be the only pin- 
ning structure in superconducting materials. Col- 
umnar disorder may also be introduced into materi- 
als for the purpose of increasing the critical current 
|13j . Far above the depinning threshold a 'moving 
Bose glass' is formed in the presence of columnar 
pins ^ . Owing to the correlated nature of the dis- 
order along the length of each vortex, the structure 
function tends to resemble that of the moving trans- 
verse glass. Furthermore, akin to the equilibrium 
Bose glass (i.e., the disorder-dominated amorphous 
structure formed by flux lines localized by corre- 
lated defects [15), the moving Bose glass displays 
a diverging tilt modulus j9|15j. Whether a unique 
voltage noise signal exists for this phase (as well 
as for the moving transverse glass) is unclear, and 
subject to this present investigation. 

The washboard frequency has been observed in 
a number of experiments [16|17|18] . However, direct 
observation of the washboard noise was achieved by 
Togawa et al. [19 , who obtained voltage noise spec- 
tra of BSCCO crystals in the mixed state subject to 
a constant current for various applied magnetic field 
strengths. For low magnetic flelds broadband noise 
was observed and attributed to plastic vortex flow. 
As the applied magnetic fleld was increased, the 
broadband noise signal reached a maximum value 
and then decreased again, while a narrowband noise 
peak that corresponded to the washboard frequency 
emerged. Upon increasing the magnetic fleld fur- 
ther, the characteristic frequency of the narrowband 
noise grew owing to a tighter flux line packing of 
the vortices, and hence a shorter vortex lattice con- 
stant. The narrowband noise signal also decreased 
in height and increased in width. The reason for 
this is apparently still not fully understood. 

Washboard noise has also been detected unam- 
biguously in a number of two-dimensional numeri- 
cal simulations; as characteristic examples, we men- 
tion the following: Olson et al. \^0f performed mole- 
cular dynamics simulations of vortices that were 
driven through a system of randomly placed de- 
fects. Upon varying the drive strength they noticed 
that the number of regimes available to the system 
above the plastic flow phase depended on the vor- 
tex interaction strength. For systems with interme- 
diate to high interaction strength and large applied 
drive, the vortices entered a coupled channel regime 
indicated by sixfold coordination of the structure. 
In this regime washboard noise was observed. Fur- 
thermore, it was noted that the washboard signal 
intensity decreased as the system size was increased. 
It is believed that this was due to multiple domains 
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forming in the vortex lattice resulting in decoher- 
ence of the noise signal. 

Kolton et al. ^21j also performed two-dimension- 
al molecular dynamics simulations to investigate 
Fiory steps [1^ for different vortex velocities. In 
order to understand the relationship between these 
steps and the temporal order in the different regi- 
mes, the voltage power spectrum was investigated 
without an applied ac drive. The authors observed 
the evolution of the power spectrum from broad- 
band noise to the emergence of a narrowband peak 
as the applied dc drive was increased, and also found 
higher harmonics. Two-dimensional computer sim- 
ulations employing Langevin dynamics and varying 
drive and pinning strength have been performed as 
well |22| . with similar washboard noise results. 

In contrast, we are aware of only one three- 
dimensional investigation of the washboard noise. 
Using an anisotropic XY model Chen and Hu TT] 
investigated the first-order transition from the mov- 
ing Bragg glass to the moving smectic. In the mov- 
ing Bragg glass phase narrowband noise that corre- 
sponded to the washboard frequency was observed 
along with harmonics for various driving strengths 
at zero temperature. The washboard peak itself per- 
sisted for T > 0, while its harmonics disappeared. 
The authors also found that the moving Bragg glass 
turned into a moving liquid at high drive, and hence 
that the washboard signal was destroyed. They ar- 
gue this to be due to thermally activated vortex 
loops inducing dislocations in the Bragg glass [11] . 

The purpose of this work 23J is threefold. First, 
we wish to establish confidence in our novel simula- 
tion approach to modeling vortex motion by quali- 
tatively comparing our results to well-known super- 
conducting vortex behavior. Second, we wish to in- 
vestigate the evolution of the narrowband noise as- 
sociated with the washboard frequency for increas- 
ing vortex density (i.e., increasing magnetic flux 
density) in the presence of randomly distributed 
point as well as columnar defects. Finally, we are 
interested in studying the effects of different types 
of pinning centers on the narrowband voltage noise. 
Below the critical current, in the presence of colum- 
nar defects, it has been predicted that vortices hop 
between pinning sites temporarily trading elastic 
deformation energy in the form of double-kinks and 
half-loops for a lower overall energy configuration 
[14] . Well above the depinning current remnants of 
these half-loop excitations and double kinks may 
still exist. At high driving values these remnant ex- 
citations would occur predominantly in the direc- 
tion of the drive, while other excitations would be 
suppressed by the localizing effect of the columnar 
pins. Obviously, the impact of such vortex excita- 
tions on the power spectrum cannot be addressed 
by a two-dimensional simulation, but require a full 
three-dimensional model. 



Based on the effective free energy for interact- 
ing magnetic flux lines in the London approxima- 
tion, and subject to attractive pinning centers [2], 
we have developed a three-dimensional Monte Carlo 
simulation code [24] to study the effects of disorder 
on the velocity / voltage power spectrum and the 
two-dimensional static vortex structure factor in 
the plane transverse to the magnetic field for driven 
vortices in the non-equilibrium steady state. Specif- 
ically, we compare results for point and columnar 
defects. We also measure the average radius of gy- 
ration in order to examine the effects of the differ- 
ent defect types on the shape or thermal 'wander- 
ing' of the elastic flux lines along the magnetic field 
(z) direction. The simulation results reported here 
should be contrasted with our earlier findings for 
non-interacting flux lines in the presence of various 
disorder distributions [Mj . 

Our results display many similar features for 
both defect types. As the vortex density is increased 
for systems with either weak point or correlated 
disorder, positional ordering is observed to increase 
in the structure factor plot. For columnar defects 
the vortex structure factor is found to change from 
that characteristic of a typical liquid, to a smectic, 
and eventually an ordered triangular lattice. For the 
case of point defects, we only observe the triangu- 
lar array in the parameter region studied here. We 
find that the structure factor plots at low vortex 
densities in the presence of point disorder appear 
qualitatively similar to the results for columnar de- 
fects at higher densities. As the structure factor be- 
gins to display positional order in the direction of 
the drive a narrowband noise signal in the velocity 
noise power spectrum corresponding to the wash- 
board effect is detected for both defect types. Asso- 
ciated with the washboard peak are harmonics, the 
ratios of which initially appear to be related to the 
type of pinning defect in the system. We present 
results that suggest these harmonic ratios are in 
fact not dependent on the type of pinning centers 
present in the sample. 

To further examine these potentially distinct ef- 
fects on the velocity or voltage power spectrum, we 
vary the effectiveness of the pinning centers. For a 
fixed vortex density the positional arrangement of 
the pinning sites in the system is changed from ran- 
domly distributed (i.e., point-like defects) to corre- 
lated along the z axis (i.e., columnar defects). The 
velocity fluctuation power spectra as well as the 
structure factor plots are examined for this series 
of simulations as well. We then compare our find- 
ings to the results obtained with increasing point 
defect pinning strength. We find that whereas the 
power spectrum and structure factor evolve in a 
qualitatively similar manner when varying the 'pin- 
ning effectiveness' through either method, there ap- 
pear marked differences in the behavior of the mean 
radius of gyration. Namely, as the point pinning 
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strength increases, vortices tend to stretch and de- 
form following a Vg (X behavior, where U repre- 
sents the pinning potential depth. In contrast, we 
observe the radius of gyration to saturate for in- 
creasing columnar defect length. The behavior is 
best described phenomenologically as Vg cx e~'°/' 
where I denotes the length of the columnar defects 
(and Iq gives the length scale in z direction). 

In the following section 2, we describe our model 
and the simulation algorithm in detail. Section 3 
contains our simulation results, as already summa- 
rized above. Finally, in section 4, we conclude and 
provide an outlook for further investigations. 



2 Model Description and Simulation 
Algorithm 

In our Monte Carlo simulations, vortices are consid- 
ered in the London approximation (with the Lon- 
don penetration depth large compared to the co- 
herence length, A » ^) as discretized elastic lines 
[l4] (see also Refs. |25|26) ). The elastic energy asso- 
ciated with the line tension of Ny flux lines is taken 
to be 



/ dz 


dviiz) 


Jo 


dz 



(2.1) 



where ri{z) describes the configuration of the ith 
vortex by specifying its two-dimensional position 
r as function of the coordinate z {0 < z < L) 
along the magnetic field direction. The line stiffness 
is given by ei = eo In(^a6/Aa6)^"^, where Xab is 
the in-plane London penetration depth, and ^ab the 
in-plane superconducting coherence length, ep — 

(0o/47rAah) sets the overall energy scale, and (j)o — 
hc/2e is the magnetic flux quantum. The expression 

p.ip for the elastic energy holds if \dri{z) / dz\ <ti 
r~^, where F'^ = Mz/M±_ denotes the effective 
mass ratio for the elastic line. In this study we 
model high-Tc materials for which ^ 1. In the 
simulation each flux line is represented by Np points 
located at (r^, Zi). Each point is confined to a con- 
stant Zi (a separate ab plane) and interacts with 
its nearest neighbors above and below via a simple 
harmonic potential. 

The total interaction energy between all pairs of 
distinct vortices is 




with the pair potential V{r) = 2eo-ft'o(r/Aaf,). Here, 
Kq is the modified Bessel function of zeroth order, 
and can be described qualitatively as diverging log- 
arithmically as r — > and decreasing exponentially 
for long distances r ^ Xab- Interactions between 



vortices occur only within the ab planes (i.e., con- 
sistent with the London limit we neglect any cross- 
plane interactions); this approximation is valid as 
long as the requirements for Eq. (|2.ip are satisfied. 
For the simulation, we consider a system of exten- 
sion Lx and Ly in the x and y directions with pe- 
riodic boundary conditions. For each vortex pair, 
we compute its minimal distance within this rect- 
angle and its periodic images adjacent to it, and 
use that distance to evaluate the interaction po- 
tential. As a consequence of this nearest-image ap- 
proximation, the interaction potential is cut off at 
distance m\Ti{Lx/2,Ly/'2,). To minimize the effects 
of the cut-off, Xab has been decreased to prevent 
numerical artifacts observed in the simulation [27j . 
We have also run Monte Carlo simulations for a 
cutoff length twice the original length used in the 
bulk of this study. To accommodate the increase in 
length the system area was increased by a factor of 
four. For the larger system the washboard noise is 
recovered for both columnar and point defects, and 
higher harmonics are observed for point defects. 

Material defects in the system are represented 
by a distribution of cylindrical potential wells. 




with Vn{r,{z)) = Y.k=i Ue{bo - \r,{z) - ri^^l). 
Here 6o is the pin radius in the ab plane, & denotes 
the Heaviside step function, r^^"^ indicates the loca- 
tion of the fcth pinning center, and U characterizes 
the well depth. 

Finally, in the presence of an external current 
the flux lines experience a force per unit length 
/l = (f'oz X J/c, therefore a corresponding work 
term is introduced: 

W^-J2 fL■r^{z)dz. (2.4) 

This work contribution favors vortex motion along 
the direction of the force while suppressing motion 
against it. For this investigation, the force is always 
applied in the x direction. The total energy of the 
system reads i?tot = El + i?int + Ed + W. 

In our simulations the applied magnetic field is 
taken parallel to the c axis (oriented parallel to z); 
therefore, at time t = straight lines are placed 
vertically in a system of size L^x LyX with peri- 
odic boundary conditions in all directions. We have 
found that the initial configurations of the individ- 
ual lines do not affect the steady state attained at 
long times. Defect centers at positions r^,^'' are also 
distributed throughout the system, either randomly 
or aligned parallel to the c axis to model columnar 
pins. The state of the system is then updated ac- 
cording to standard Monte Carlo Metropolis rates. 



T.J. BuUard et al.: Vortex Washboard Voltage Noise in Type-II Superconductors 



5 



When the number of attempted updates is equal to 
the number of points that make up a flux line mul- 
tiplied by the total number of lines in the system, 
this constitutes a single Monte Carlo step (MCS) 
and serves as the unit of time in the simulation. 
While this 'driven diffusive system' simulation ap- 
proach, introduced by Katz, Lebowitz, and Spohn 
[28] , is likely best suited to model thermally acti- 
vated motion close to equilibrium, we find that our 
results in the driven regime, considerably away from 
equilibrium, are quite similar to those observed in 
other studies, as mentioned above. 

Nevertheless, it is not at all obvious even for 
steady states far from thermal equilibrium which 
choices of 'microscopic' Monte Carlo rates yield the 
most 'realistic' representations of an experimental 
system, even in a coarse-grained view. Gotcheva 
et al. |10] have recently brought into question the 
general validity of the driven diffusive Metropolis 
Monte Carlo dynamical method in their simulations 
for vortices on a two-dimensional lattice. The non- 
equilibrium steady states obtained using Metropo- 
lis Monte Carlo and continuous-time Monte Carlo 
dynamical rules, respectively, were compared as a 
function of temperature and driving force. The re- 
sults differed dramatically depending on which dy- 
namical rule was used. In some instances at least, 
the Metropolis algorithm yielded a spatially disor- 
dered moving steady state while the continuous- 
time Monte Carlo rules preserved positional order 
(in finite-size systems) over much of the non-equili- 
brium phase diagram. The authors argue convinc- 
ingly that the lack of order observed in the Metro- 
polis simulation was due to intrinsic randomness in 
the updating rules. It would certainly be interest- 
ing and worthwhile to probe to what extent these 
findings also apply to three-dimensional off-lattice 
simulations, which are presumably less likely to be 
stuck in long-lived metastable configurations. In our 
current study we find that spatial order survives in 
the non-equilibrium steady state for an extended 
range of flux density values. 

Quite generally, it is crucial for the analysis of 
out-of-equilibrium systems to carefully investigate 
alternative approaches to the description of their 
dynamics in order to probe their actual physical 
properties rather than spurious artifacts inherent in 
any mathematical modeling. Different mathemati- 
cal and numerical representations of non-equilibri- 
um systems in fact rely on various underlying a- 
priori assumptions that must be tested a-posteriori 
by comparing their respective results. For example, 
many computer simulation studies of driven flux 
line systems invoke stochastic Langevin equations, 
wherein one assumes that all fast degrees of freedom 
are aptly captured in terms of uncorrelated white 
noise (see, e.g., Refs. j20|29|22j ). Such a mesoscopic 
representation of the dynamics is usually adequate 
in thermal equilibrium where the form of the input 



noise correlations is severely restricted by fluctua- 
tion-dissipation relations. Yet the large-scale and 
long-time properties of nonlinear Langevin stochas- 
tic differential equations away from thermal equilib- 
rium are well-known to be often drastically affected 
by the functional form or even strength of the as- 
sumed noise correlations, which are not uniquely 
determined by Einstein relations any more. One 
would, e.g., suspect that the ffux line correlations 
along the magnetic field axis should be reflected in 
the noise spectrum and relaxational features. It is 
therefore imperative to test a variety of different nu- 
merical methods and compare the ensuing results in 
order to identify those properties that are generic 
to the physical system under investigation. 

The parameter values used in our simulations 
correspond to typical high- Tc materials, and are re- 
ported in units of 60 (the pin radius) and interaction 
energy scale eg. The parameters ^ab, ei, and are 
chosen to be 0.5 6oj 0.25 eo, and 16, respectively. In 
this study we examine vortex motion in the weak 
pinning regime. The pinning potential U has been 
given a value of Uq ~ 0.03125 eo while the predicted 
value is approximately 0.5 eg JJJ. The penetration 
depth Xab is assigned a value of 16 bo which is about 
1/3 of typical high- Tc values. As previously men- 
tioned this choice was made in order to minimize 
artifacts due to the interaction cut-off. The average 
separation distance between randomly distributed 
defects in each ab plane is taken to be 15 6o- The 
maximum distance for a point on any flux line to 
move is limited to 6o/2, to help avoid vortices 'hop- 
ping' over defects. We have used a discretization 
along the fleld direction of = 20 parallel planes. 

As is well-known in flnite-size vortex simula- 
tions, a square planar geometry favors ordering into 
a square lattice whose conflgurational energy is only 
slightly above that of the triangular Abrikosov lat- 
tice pO]. We have thus chosen the system's planar 
aspect ratio L,^ : Ly = 2 : -\/3 in order to eas- 
ily accommodate a triangular lattice. The vortices 
are then placed in the system prearranged in a tri- 
angular array with a lattice vector oriented along 
the x-direction. Choosing this aspect ratio allows 
an even square number of vortices to 'flt' in the 
system while arranged in this configuration. In the 
simulation runs, the vortex lattice arranged in lo- 
cal energy minima configurations by either align- 
ing a lattice vector parallel to the system's hori- 
zontal axis, rotating the lattice vector by 30° from 
the horizontal, or by arranging such that the lattice 
twisted about the system at a chiral angle satisfying 
periodic boundary conditions. When in our simula- 
tions an external force is applied and the vortex 
system driven, we find that the lattice maintains 
its initial orientation. However, by approximately 
doubling the defect pinning strength the lattice re- 
orients such that the principal lattice vector points 
in the direction of the applied drive. Drive-induced 
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reorientation has been observed in experiments [31j 
and simulations [22]. 

Since the present study is primarily concerned 
with the effect of defect correlations on the dynam- 
ics, the temperature is chosen such that T/T* < 1, 
where T* = fc^^VeiCTo^o is the temperature above 
which entropic corrections due to thermal fluctu- 
ations become relevant for pinned flux lines [T3]. 
Here, thermally induced stretching and wandering 
of the flux lines are largely suppressed (as long as 
the vortices remain pinned), and the results can be 
interpreted in terms of low-temperature kinetics. 
Our simulations are thus usually run at k^T per 
unit length equal to 0.004 eq- The average or cen- 
ter of mass (CM) velocity for each vortex is then 
calculated and averaged over all vortices: 



— T 



rc-m,{t + r) - rcmi(i) 



(2.5) 



where rcmi denotes the center-of-mass position of 
the ith vortex, and r is the time interval between 
measurements, r is set to 30 MCS, and the simula- 
tion is then run for 10^ MCS to arrive at a steady 
state. Data are subsequently collected for the next 
2.5 X 10^ MCS. 

From the collected data, we obtain the two-dim- 
ensional static structure factor in the plane trans- 
verse to the magnetic field, 



^(k) - / (p(O)p(r)) e 



-ikr 



dr , 



(2.6) 



where (p(O)p(r)) denotes the density-density corre- 
lation function, for the driven vortices in the non- 
equilibrium steady state. We also measure the av- 
erage radius of gyration. 



1 



1/2 



(2.7) 



in order to examine the effects of the different defect 
types on the shape or thermal 'wandering' of the 
elastic flux lines along the magnetic field (z) direc- 
tion. Additional information about the detailed dy- 
namics is encoded in the velocity fluctuation power 
spectra 



v{t) e'"'* dt 



(2.8) 



These velocity power spectra have been appropri- 
ately windowed in order to minimize spectral leak- 
age. Since vortex motion across the sample induces 
a voltage drop, the velocity fluctuations are exper- 
imentally directly accessible as the measured volt- 
age noise. Each of the above observables is mea- 
sured for various vortex densities and defect config- 
urations. Vortex densities are reported as a count 
of the number of vortices in a unit system of area 



A = 150^ bo X 150 6o- The vortex structure factor 
and velocity noise plots are averaged over time and 
typically over 40-50 disorder realizations. 



3 Monte Carlo Simulation Results 

3.1 Current-Voltage (I-V) Characteristics 

In order to validate our model and simulation al- 
gorithm we have examined the average vortex ve- 
locity as a function of the applied force and com- 
pared our results to well-established experimental 
and numerical findings. In experiments, the driv- 
ing force is proportional to the externally applied 
current, and the induced voltage across the sample 
proportional to the mean flux line velocity, hence 
the current-voltage (TV) characteristics are given 
by the force-velocity curves in our simulations. 

Our simulation results for systems with random- 
ly distributed columnar pins and point defects are 
plotted for various flux densities in Fig. 13.11 In ad- 
dition to averaging over all flux lines in the sys- 
tem, the data are averaged over time (25,000 Monte 
Carlo steps) and five disorder realizations. For both 
graphs the number of effective pinning sites is equal, 
as are the pinning strengths (per unit length) Uq for 
all the defects. 

The vortices are seen to remain pinned for both 
columnar and point defects up to a critical depin- 
ning force, beyond which the system takes on an 
average velocity that increases with increasing ap- 
plied drive. As expected the depinning threshold is 
higher for columnar defects [Fig. 



3.1(a) since the 



attractive force from each pinning site adds coher- 
ently over the length of the flux line. The inset in 
Fig. 3.1(b) depicts the region of the curve associ- 



ated wltlTdepinning for point defects, showing just 
the data for the densest and most dilute systems. 
We observe (again as expected) that denser vor- 
tices depin at lower driving currents for either type 
of pinning center; estimates for the depinning cur- 
rent for columnar d efects a s obtained from our data 
are depicted in Fig. 3.1(c) As the repulsive interac- 



tions between neighboring vortices grow with flux 
density, this helps the vortices to overcome the de- 
fect pinning potentials. We find that these results 
are qualitatively comparable to experimental find- 
ings in high-temperature superconducting materi- 
als (see, e.g., Refs. [32"33"34|), and in accord with 
results obtained by means of Langevin molecular 
dynamics simulations |35| (in the three-dimensional 
regime, see also Ref. [35]), giving us confidence in 
the validity of our model and algorithm. 

However, we note two artifacts in the TV re- 
sults occuring at higher driving current and large 
vortex densities. First, while at large currents linear 
Ohmic behavior should ensue, we find that the ve- 
locity values saturate at high applied drive values. 
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Fig. 3.1. Velocity (induced voltage) vs. force (applied 
current) curves (I-V characteristics) for (a) randomly 
distributed columnar pins, and (b) point defects. For 
both defect types the depinning threshold decreases 
with increasing vortex / flux density. This can be seen 
most clearly for point defects in the inset in (b), which 
amplifies the low current regime: here, the force values 
on the X axis of the inset range from to 0.0025, the 
velocity ranges from to 0.0016 and is given in units of 
b(i/ MCS (pin radius per Monte Carlo step). The follow- 
ing symbols represent vortex densities reported as the 
areal density p in a fixed area A — 150-^ feo x 150 6o: 
A - 144/y4 = 0.00554/6g, ■ - lOO/A = 0.00385/feg, 
* - QA/A = 0.00246/6o, O - 36/A = 0.00139/6g, 
X - IG/A = 0.00062/6o- (c) Estimate of the depinning 
force fc as function of the vortex density p for random 
columnar defects. The inset depicts how the I-V curve 
for p = 14:4/ A saturates at high drive values. In all plots, 
the data points are connected as a guide for the eye. 



see inset of Fig. 3.1(c) This is due to a maximal 
step size limitation imposed to avoid vortices 'hop- 
ping' over pinning defects in the system. Second, 
the TV curves are observed to cross at larger den- 
sities due to the fact that the system is updated 
locally rather than globally, in the Metropolis algo- 
rithm. At higher vortex densities larger local moves 
are suppressed by the repulsive potential of nearest 
neighbors resulting in a lower average velocity and 
an I-V curve with a lower slope than that of lower 
density systems. 

In the following, we mostly report data that 
have been collected at a driving force / = 0.04. 
This high value was chosen so that washboard noise 
could be observed over a range of accessible vortex 
densities; however, in some instances the TV curve 
has already begun to saturate. Since we are study- 
ing the velocity noise spectrum relative to the mean 
velocity, we do not expect specific artifacts caused 
by the saturation in the voltage noise, nor do we ob- 
serve any as compared to results obtained at lower 
drive values / = 0.01 (sec the discussion at the end 
of SccEXT]). 

3.2 Narrowband Noise Characteristics 

Power spectral density plots and intensity plots of 
the structure factor 5'(k) have been obtained for 
several vortex densities as shown in Figs. 13.31 and 
13.81 Unless otherwise indicated, all power spectral 
density plots are double-logarithmic with the fre- 
quency on the X axis ranging from 0.001 to 0.1 rad 
/ MCS and power in normalized units on the the 
y axis in the range 2.5 x 10"^" . . . 1 x 10~*. Dis- 
tinct peaks are observed as well as higher harmon- 
ics in many of the power spectra for columnar de- 
fects. Peaks are observed in all of the power spec- 
tral plots for point defects. In both cases harmonics 
are always located at integer multiples of the fun- 
damental frequency. In Fig. 13.21 the fundamental 
frequencies are plotted versus the number of vor- 
tices per unit system size along with the predicted 
washboard frequencies for both point and columnar 
disorder. The washboard frequency is calculated by 
simply dividing the measured average vortex veloc- 
ity {vx) (parallel to the applied drive, and averaged 
over time and defect realizations) by the vortex tri- 
angular lattice constant. Due to the aspect ratio 
of the system this distance is obtained by dividing 
the system length in the direction of the drive {Lx) 
by the square root of the number of vortices in the 
system; hence, 

Lo = 2T:{vx)y^/Lx. (3.1) 

We obtain good agreement between the measured 
and predicted values indicating that the fundamen- 
tal frequency in the power spectra plots is indeed 
the washboard frequency. Error bars for the mea- 
sured frequencies are estimated by taking the full 
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Fig. 3.2. Washboard frequency versus vortex density p 
for (a) point and (b) columnar pinning centers. Vortex 
density is reported as a count of the number of flux 
lines in a fixed area of size A — 150-^ 6o x 150 6o. The 
results show good agreement between the measured and 
calculated values for both defect types. The frequency 
uo is given in units of rad / MCS. Error bars for the 
measured values are obtained from the full width at 
half- maximum of the washboard peaks in Figs. 13.31 and 
13.81 Error bars for the calculated values are smaller than 
the data points. ▲ - measured, ■ - calculated. 



widths at half maximum of the fundamental peaks, 
while the uncertainties for the calculated values are 
obtained from the standard deviations of the aver- 
age velocities. 



3.2.1 Randomly Distributed Columnar Defects 

Noise characteristics and structure. Results for flux 
lines of different densities interacting with random 
columnar defects are displayed in Fig. 13.31 With 
the average spacing between defects set to 15 6o the 
number of columns in a unit area A = 150-^ 6o ^ 

150 6o is 115. For a flux density p = 16/ A (a filling 
fraction of ~ 1/7) only broadband noise is obse rved 
in the velocity power spectrum, see Fig. 



3.3(a) The 



corresponding diffraction pattern shows a ring indi- 
cating an isotropic liquid or amorphous solid. A typ- 
ical Delaunay triangulation for the vortex positions 
in a snapshot of a particular run with 16 lines is de- 
picted in Fig. 3.4(a) The plot shows the presence 
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Fig. 3.3. Velocity / voltage power spectra measured 
in the direction of the drive {x direction) and struc- 
ture factor plots for increasing vortex density in the 
presence of randomly distributed columnar defects, (a) 
For a density p = 16/^4 — 0.00062/6o an isotropic liq- 
uid is observed along with a broadband noise signal in 
the power spectrum, (b) At p = 36/A = 0.00139/bg 
two (off-center) peaks appear in the structure factor lo- 
cated transverse to the drive direction, accompanied by 
a narrowband signal in the velocity power spectrum, 
(c) Six peaks appear in the structure factor plot for 
p = 64/A = Q.00246/6g. The intensities of the peaks 
with a wave vector component in the x direction are 
lower than those with only a y component. In the veloc- 
ity power spectrum the washboard frequency peak nar- 
rows indicating greater temporal coherence, and higher 
harmonics become more visible, (d) As the density is 
increased to p = 100/A — 0.00385/6o the structure be- 
comes more fully ordered into a triangular array, and 
the washboard peak sharpens. 



of a number of topological defects in t he vorte x sys- 
tem. The 'time exposure' plot in Fig. 3.4(c) of the 
ffux line center-of-mass positions produces trajecto- 
ries reminiscent of the 'braided rivers' observed in 
two-dimensional plastic flow simulations |29j . The 
trajectories appear to move and cross within wind- 
ing channels. For a density p = 36/ A the two peaks 
located perpendicular to the drive direction fully 
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(a) p = 1&/A 



(b) p = 64/yl 




(c) p = 1&/A 



(d) p = 64/A 



Fig. 3.4. Delaunay triangulation plots for the positions 
of (a) 16 and (b) 64 flux lines per area A in the pres- 
ence of random columnar defects. Topological defects 
are marked. A disordered structure is obtained for flux 
density p = 16 /A, while a triangular array free of topo- 
logical defects is found for p = 64/ A. A 'time exposure' 
for densities 16/^ and 64/^ is plotted in (c) and (d) 
to illustrate the flux line motion. Open circles indicate 
columnar defect locations, while black lines trace the 
trajectories of the average position of each vortex. For 
p = 16 /A, the vortex trajectories cross suggestive of 
plastic flow, while for p = 64/A parallel channels form. 



emerge, suggestive o f the p redicted moving trans- 
verse glass, see Fig. 



3.3(b) Small intensity peaks 



with wave vector components paraUel to the drive 
are also visible. 

As the vortex density p is increased to 64/ A, 
peaks with an x component in the diffraction plot 
emerge more prominently and sharpen as shown 
in Fig. 3.3(c) Just as for p — 36/A, these peaks 
are smaller than those located perpendicular to the 
direction of the drive. I n a typical associated De- 
launay plot, Fig. 3.4(b) topological defects in the 



vortex lattice disappear, and the flux line trajec- 
tories form parallel channels, running in compara- 
tively straight lines, Fig. |3.4(d)[ At p = 100 /A the 
structure plot reveals a well-ordered lattice of flux 



hues, see Fig. 3.3(d) 



These results are a good demonstration of the 
competing energy scales in the system. The spa- 
tially randomly distributed pinning sites favor dis- 
ordered flux line structures while the vortex repul- 
sion induces ordering into a regular array. For p = 
16 /A the structure factor displays a random vor- 
tex configuration, indicating that the lattice struc- 
ture is dominated by the columnar pinning sites; 
indeed, the system is close to the depinning thresh- 
old, see Fig. 



3.1(c) 



Individual vortices are pinned 
for periods of time that are long compared to the 
time it would take for the flux lattice to move one 
lattice constant, thus preventing the formation of 



1000 
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Fig. 3.5. Washboard peak power for increasing vor- 
tex density for columnar (square data points) and point 
defects (triangles). In either case the peak intensity is 
observed to decrease as the flux density is increased. 



the regular structure. At density 36/A, the system 
has moved away from the depinning threshold, and 
the repulsive forces begin to separate the vortices 
into parallel channels resulting in spatial periodic- 
ity in the y direction (perpendicular to the drive), 
and by reaching p — 64/ A the system of moving 
vortices is dominated by the repulsive vortex inter- 
action potential. At this stage individual channels 
begin to couple as additional periodicity appears in 
the X (drive) direction. Distinct peaks emerge in the 
diffraction plot with the peaks possessing a wave 
vector component in the x direction smaller than 
those with only a y component due to weaker cou- 
pling. As the density is increased to p = 100/ A the 
repulsive energy between the vortices grows further 
leading to a stronger coupling between transverse 
channels and eventually a symmetric triangular lat- 
tice structure emerges. 

The spatial structure of the vortex array is re- 
flected in the associated velocity or voltage power 
spectrum. At p = 16/A a broadband signal is visi- 
ble corresponding to the isotropic vortex distribu- 
tion and their consequent random incoherent mo- 
tion shown in Fig. 3.4(c) Similar broadband noise 
that is associated with incoherent flux transport has 
been experimentally and numerically observed in a 
number of studies [29 6 19J. For p = 36/A, as the 
vortices begin to arrange into a lattice a narrow- 
band peak appears corresponding to the x wave vec- 
tor components that appear in the structure factor 
plot. At a density 64 /A the washboard signal nar- 
rows as the vortex system becomes more structured. 
A sharper peak indicates a greater temporal coher- 
ence between the vortices which is not surprising; as 
the density of the vortices is increased each vortex 
is more tightly held in its place and each line 'stiff- 
ened' by the repulsive forces exerted by its neigh- 
bors. For similar reasons the remnants of the broad- 
band noise ffatten out more. This trend continues 
as the flux density is further increased to 100 /A. As 
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the density is increased from p = 6A/A upward we 
observe a decrease in the washboard peak power, as 
shown in Fig. 13.51 The decrease of the fundamental 
peak intensity is due to the increase in stiffness of 
the lattice structure with growing vortex density. 
As each flux line passes through a defect site it will 
be less affected by the pin in a denser system re- 
sulting in smaller velocity fluctuations and hence 
reduced power output. 

Washboard peak harmonics. In Table [XT] the ra- 
tio of the intensities of the first and second har- 
monic with respect to the fundamental peaks are 
recorded for the velocity noise measured in the x 
and y directions. These ratios vary as the flux den- 
sity increases indicating a change in the shape of the 
velocity vs. time trace. It is observed that the ra- 
tios approach values similar to those measured for 
point defects (listed in Table [3721 below). We will 
investigate the relationship between the underlying 
defect type and the harmonics ratios in Sec. 13.3.11 
A narrowband signal is also measured in the y di- 
rection transverse to the drive for densities of 6A/A 
through 144/A flux lines (not shown). The ratios 
of the harmonics to the fundamental peak intensity 
in the y direction are also reported in Table 13.11 
The frequency of the fundamental peak as well as 
the higher harmonics are identical to the frequen- 
cies measured in the x direction. While the power is 
lower in the y direction (by factors of 10 to 40), the 
ratios of the peak intensities are similar. Since there 
is no drive in the y direction, these observations 
could perhaps be interpreted as a suppression of 
transverse motion occurring at the same frequency 
as the washboard motion. One likely explanation is 
that as a vortex becomes temporarily trapped in 
a pinning potential, fluctuations transverse to the 
motion are suppressed until the vortex departs from 
the pinning site, resulting in periodic behavior. 

Radius of gyration. In order to obtain additional 
information about the three-dimensional shape of 
the elastic vortex lines moving through a sample 
with columnar defects, we have obtained radius of 
gyration data, averaged over vortices, time, and de- 
fect configurations. We interpret the radius of gyra- 
tion as the distance the elastic line is stretched from 
its average (center-of-mass) position. The compo- 
nents of the radius of gyration in the x and y di- 



rection vs. fiux density are plotted in Fig. 3.6(a) 
As the flux density is increased the radius of gyra- 
tion decreases indicating that the vortex lines are 
straightened at higher densities due to the stronger 
repulsion exerted by their nearest neighbors. The 
data also display anisotropy in the flux line 'stretch- 
ing'. The magnitude of the radius of gyration is 
markedly greater along the direction of the drive 
than perpendicular to it, with the x and y com- 
ponents approaching each other only at quite large 
densities. For dilute systems we would expect the 
competition between the drive and the pinning po- 
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Fig. 3.6. Components of the mean radius of gyration 
for (a) columnar and (b) point pins. For both defect 
types the radius of gyration decreases with increas- 
ing vortex density. At higher flux densities vortices are 
caged by nearest neighbors, suppressing elastic flux line 
stretching. Lengths are reported in units of the pin ra- 
dius bo. k - X component, ■ - y component. 



tential to result in flexible vortices depinning in 
sections, with some parts of the flux lines leaving 
the columnar pins via 'double-kink' and 'half-loop' 
saddle-point configurations [T?'. In the presence of 
the drive the free flux line sections move forward 
and grow, stretching the vortex until it depins com- 
pletely. Since there is no drive in the y direction, the 
y component of the radius of gyration is smaller; 
however, the same trend is observed with increasing 
density as with the x component. At high densities 
the vortices are stiffer due to a smaller vortex sepa- 
ration distance. Rather than stretching, depinning 
tends toward an 'all or nothing' process with the 
tightly packed vortex ensembles effectively becom- 
ing two-dimensional. 

Comparison with results obtained at small drive. 
The power spectra results reported above have been 
obtained at a rather high applied drive / = 0.04 
chosen so that the washboard peak could be ob- 
served over a range of vortex densities. However, 
as discussed above, the velocity in the I-V charac- 
teristics (Fig. 13. ip begins to saturate at such high 
drive values, especially for systems with higher flux 
density. We believe that this is an artifact of the 
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Table 3.1. Ratio of the three largest peaks observed in the vortex velocity power spectrum in the presence of 
randomly distributed columnar defects. For each vortex density the ratios of the intensities of the second and 
third peak to the first are reported, for measurements taken in both the x and y directions. The number of runs 
over which the power spectral density plots were averaged is also listed. The corresponding power spectra peak 
ratios for a periodic piecewise linear (sawtooth) signal is included for comparison. 



vortex density p 


runs 


ratio {x direction) 


ratio (y direction) 


64:/ A = 0.00246/6o 


115 


1 : 0.20 ±0.02 : 0.030 ± 0.003 


1 : 0.25 ± 0.03 : 0.079 ± 0.009 


100/A = 0.00385/6;^ 


44 


1 : 0.31 ±0.04 : 0.055 ± 0.008 


1 : 0.30 ± 0.05 : 0.070 ± 0.010 


lU/A = 0.00554/6^ 


44 


1 : 0.38 ±0.08 : 0.120 ±0.030 


1 : 0.36 ±0.06 : 0.055 ± 0.008 


f{x) = X, < X < 27V 




1 : 0.25 : 0.11 





limited step size in the simulation. In order to in- 
vestigate whether the velocity saturation influences 
the velocity fluctuations we now examine the power 
spectra at a lower drive value / — 0.01 and compare 
the results to those for / = 0.04. We consider here 
a dense vortex system with p — IM/A subject to 
randomly placed columnar defects. This flux den- 
sity was chosen to preserve the washboard peak at a 
lower drive. Results are displayed in Fig. 13.71 These 
data suggest that the velocity saturation does not 
adversely affect the velocity power spectra. The ex- 
pected drop in washboard frequency corresponding 
to the decrease in applied drive is observed. The full 
peak widths at half maximum of the fundamental 
are measured at high and low drive and the results 
are found to be comparable (0.00036 rad/MCS* and 
0.00041 lad/MCS, respectively) indicating that the 
temporal correlations are not affected by the satura- 
tion. Both spectra display a similar harmonic ratio 
1 : 0.40 ± 0.05 : 0.10 ± 0.01. Likewise the struc- 
ture factors are found to be similar. These results 
suggest that data obtained at high drive values do 
not possess any pronounced artifacts due to velocity 
saturation. 
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Fig. 3.7. Voltage power spectra for (a) low (/ = 0.01) 
and (b) high (/ — 0.04) driving forces at p = 144/A in 
the presence of columnar defects. 



3.2.2 Randomly Placed Point Defects 

The results for the velocity power spectra and asso- 
ciated structure factors for vortices driven through 
randomly distributed point defects show many qua- 
litative similarities to those obtained for columnar 
pins, see Fig. 13.81 However, the structure factor plot 
associated with an isotropic liquid is never observed 
in our simulations, even at the lowest flux densi- 
ties. The intensities of the structure factor peaks 
increase with growing flux density, indicating an 
increasing degree of positional order in the system. 
For the lowest density system. Fig. 3.8(a) the struc- 



ture factor plot implies greater spatial periodicity 
perpendicular to the direction of the drive than 
parallel to it, similar to the case of columnar de- 
fects, suggesting that parallel channels are not well- 
coupled in this regime. This is further evident in the 
Delaunay triangulation plot shown in Fig. 3.9{a] 



Compared to the system with a den sity 16 /A, the 



Delaunay plot for p = i6/A, Fig. 3.9(b) yields 
greater alignment perpendicular to the drive ow- 
ing to the disappearance of topological defects in 
the vortex lattice. The reason for this is similar to 
that for columnar defects; the pinning sites intro- 
duce shear between parallel channels resulting in 
less-ordered rows of vortices that run perpendicu- 
lar to the drive. As the density increases the vortex 
repulsion becomes the dominant energy in the sys- 
tem, and the coupling between parallel channels is 
enhanced. 

In the velocity power spectrum for flux density 
p = 16 /A depicted in Fig. 3.8(a) the narrowband 



signal rests on a low broadband background signal. 
As the vortex density increases the broadband com- 
ponent diminishes, as do the power and width of the 
washboard peak. In general, the peak widths are 
smaller than those observed for columnar defects, 
as to be expected for uncorrelated and less efficient 
pinning centers. In the spectra where peaks are re- 
solvable, the fundamental peak and higher harmon- 
ics are located at frequencies that are identical in 
the X and y directions. The ratio of the fundamen- 
tal to the harmonics in either direction are listed in 
Table l3T2l These ratios for flux densities p = 16 /A 
through 100/j4 lines are comparable suggesting a 
similar shape of the velocity-time traces. For a sys- 
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Table 3.2. Ratio of the three largest washboard peaks observed in the velocity power spectrum for point defects. 
For each vortex density the ratios of the intensities of the second and third peak to the first are reported, for 
measurements taken in both x and y directions. The number of runs over which the power spectral density plots 
were averaged is also included. 



vortex density p 


runs 


ratio (x direction) 


ratio {y direction) 


16/ A = 0.00062/b;5 


44 


1 : 0.67 ±0.09 : 0.27 ± 0.02 


n / a 


36/ A = 0.00139/60 


44 


1 : 0.55 ±0.06 : 0.22 ± 0.02 


1 : 0.4 ±0.1 : 0.40 ± 0.10 


64:/ A = 0.00246/6^ 


276 


1 : 0.58 ±0.03 : 0.24 ± 0.01 


1 : 0.51 ±0.04 : 0.16 ± 0.01 


lOO/A = 0.00385/bg 


88 


1 : 0.56 ±0.05 : 0.25 ± 0.02 


n / a 


144/A = 0.00554/6^ 


88 


1 : 0.80 ±0.10 : 0.19 ±0.02 


1 : 0.43 ± 0.06 : 0.08 ± 0.01 
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Fig. 3.8. Power spectra and structure factor plots for 
increasing vortex density p in the presence of point 
pins, with (a) p = 16/ A = 0.00062/6g, (b) 36/ A = 
0.00139/6g, (c) 6A/A = 0.00246/bi, and (d) 100/ A = 
0.00385/feo. Narrowband washboard noise peaks are ob- 
served in all thespectral density plots, and six-fold coor- 
dination in the corresponding vortex structure factors. 



tern with p ~ lAA/A the ratio is however markedly 
different. At this point it is not clear to us which 
physical mechanism dictates the detailed values of 
these peak intensity ratios. Yet, as we shall further 
argue in Sec. 13.3.1] below, we believe that these ra- 
tios characterize the overall deformation of the vor- 
tex lattice rather than reflect the geometry of the 
pinning centers. 



The mean flux line radius of gyration (averaged 
over vortices, time, and pinning site distributions) 
vs. vortex density for point defects is plotted in 
Fig. 3.6(b) The behavior for both components is 



similar to that of a system subject to randomly 
placed columnar defects. A larger radius of gyra- 
tion is observed for low-density systems, while it 
decreases for both x and y components as the flux 
density is increased. Unlike for columnar defects, 
the pinning force for uncorrelated point disorder 
does not add coherently over the length of the flux 
lines. Whereas the random spatial distribution of 
point pins promotes thermal flux line wandering, 
the stretching of the vortices while moving through 
a defect is not as severe. Comparing the data to 
those for columnar defects with identical flux den- 
sities, we note that the magnitudes of both com- 
ponents of the radius of gyration are smaller for 
point defects. As the density is increased both com- 
ponents of the radius of gyration for point defects 
approach the same values as also seen in the corre- 
sponding curves for columnar defects. 

We have also examined the effect of the total 
vortex length on the velocity power spectrum. 
Data for — 60 (three times the length of the 
vortices shown in all previous results) at a den sity 
of 64 lines have been included in Figs 



3.10(a) 



20 in norma 



and 
ized 



3.10(b) (The peak power for L 
units is lAA-10-^±9- 10"^ compared to 4.2 • lO-'^i 
3 • 10~^ for Lz ~ 60.) Qualitatively, the findings are 
quite similar to those for a shorter vortex length; 




Fig. 3.9. Delaunay triangulation for vortices with den- 
sities (a) p — 16/A and (b) 36/A, interacting with 
randomly distributed point pinning centers. Topologi- 
cal defects are indicated. 
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Fig. 3.10. Velocity fluctuation power spectra and 
structure factor plots for a flux density p = 64/A for 
(a) columnar and (b) point defects in a spatially ex- 
tended system with Lz = 60, i.e., all flux lines have 
been lengthened by a factor of 3 compared to the pre- 
vious runs. The results for columnar defects shown in 
(a) are very similar to those for the same density with 
shorter lines, as depicted in Fig. 3.3(c) On the other 
hand, for point defects (b) the power drops noticeably 
compared to Fig. 3.8(c) 



however, in the presence of point defects, compared 
to the results shown in Fig. 3.8(c)| the narrowband 



power is lower. On the other hand, a similar inten- 
sity drop is not observed in the case of columnar 
defects, compare Fig. 3.10(a) to Fig. 3.3(c) The ad- 



ditional flux line length and the lack of spatial cor- 
relation in the z direction for point disorder further 
demonstrates the difference in pinning efficiency be- 
tween point and columnar defects. For columnar 
pins the lengths of both the defect and the vortex 
span the height of the sample just as they did for 
the shorter system; hence we observe a similar effect 
on the motion and the power spectrum. For point 
defects, as the length of the vortex is increased the 
effect of a single point defect on the longer vortex 
decreases as a whole resulting in a smaller power 
output: local fluctuations are averaged out. We note 
that naturally these distinctions are only observable 
in three-dimension simulations. 



3.3 Variable Pinning Efficiency 

In this section we examine the effect of varying pin- 
ning efficiency on the vortex structure factor and 
velocity fluctuation power spectrum. We shall study 
the cases of different point pin strenghts and linear 
defects of fixed strength but varying length, thus 
interpolating between uncorrelated point disorder 
and correlated columnar defects. We are particu- 
larly interested in whether point and columnar de- 
fects display different washboard harmonic signa- 
tures in the velocity power spectrum. 




0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 
U/Un 
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Fig. 3.11. (a) Velocity spectrum washboard peak 
power for increasing point pinning strength U . The vor- 
tex density is held constant at p = 64/A. As the pinning 
strength is increased the washboard signal grows to a 
maximum value around 1.7 f/o (b) and then decreases 
as the system transitions to a broadband signal around 
1.9 f/o (c). 



3.3.1 Variable Point Pinning Strength 

To carry out this investigation simulation runs were 
performed at fixed fiux density p = 64 /A while 
the point pinning strength was increased from the 
value Uq used before up to i7 = 1.9J7o- For each 
pinning strength results were obtained by averag- 
ing over random distributions of point defects. The 
number of pinning sites per run was held constant; 
only the pinning strength was changed. The power 
of the washboar d peak is plotted versus pinning 
strength in Fig. 



3.11(a) We note that the power 



reaches a maxim um value at approximately 1.7 Uq, 
see Fig. 



3.11(b) 



The increase in power is due to 
the stronger pinning causing larger velocity fluctu- 
ations in the vortex array. The width of the peak 
increases along with its power as indicated by the 
error bars in Fig. 3.11(a) The signal then decreases 



as the power output transitions to a broadband sig- 
nal coinciding with decreasing coherent vortex mo- 
tion. Fig. 



3.11(c) 



Washboard harmonic peak ratios for a few pin- 
ning strenghts are recorded in Table 13.31 We com- 
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Table 3.3. Intensity ratio of the three largest velocity 
power spectrum peaks for increasing point defect pin- 
ning strength, recorded as multiples of the initial pin- 
ning strength Uo . For each pinning strength the ratios of 
the intensities of the first and second peak to the third 
are given for measurements taken in the drive (x) direc- 
tion. The number of runs over which the power spectral 
density plots were averaged is also included. 



U/Uo 


runs 


ratio (x direction) 


1.1 


44 


1 : 0.55 ±0.06 : 0.150 ± 0.020 


1.3 


44 


1 : 0.33 ±0.04 : 0.045 ± 0.008 


1.5 


44 


1 : 0.20 ±0.03 : 0.012 ± 0.004 



pare these results to the columnar defect data in 
Table 13.11 and note that similar harmonic ratios are 
observed for the weaker columnar defects at a den- 
sity of p = 64/ A and the stronger point defects 
with U — 1.5 C/q- These findings demonstrate that 
the harmonic ratios are not unique signatures of 
the material defect types compared in this study. 
We speculate that the ratios are rather dependent 
on the amount of deformation of the vortex lattice 
by the defects regardless of whether or not the pins 
extend along the length of the vortices. 



3.3.2 Variable Linear Defect Length 

To further compare the effects of point and colum- 
nar defects on the velocity power spectrum, we have 
investigated systems with varying columnar defect 
lengths at a constant vortex density p — 16 /A. Each 
set of results is obtained by averaging runs over ran- 
dom distributions of linear defects of a particular 
length I, with the total number of pinning sites held 
constant. Defect lengths vary from a single pinning 
site, i.e., random point defects, to a length of 10 con- 
tiguous pinning sites. (Recall that a columnar de- 
fect extending over 20 pinning sites spans the entire 
system height — 20.) The driving force remains 
at / = 0.04. As the defect length is increased the 
spatial order in the system decreases as indicated by 
the structure factor peaks in Fig. 13.121 Peaks with 
wave vector components along the drive decrease 
in intensity at shorter defect lengths compared to 
those perpendicular to the drive. With increasing 
length of the linear pins (growing defect correla- 
tions) the peak intensities become diminished. Re- 
sults for defect length Z = 1, 3, and 5 are shown 
in Fig. 13.121 By length I = 10 the washboard peak 
disappears entirely and is replaced by a broadband 
signal (not shown). We note a second peak with 
a frequency corresponding to the system size also 
appears in the power spectrum 1 — 3. While the 
cause of this peak is unclear, one possible expla- 
nation is that it originates in some type of persis- 
tent deformation in the vortex lattice. With peri- 
odic boundary conditions this lattice deformation 




0.001 0.01 0-1 

VJ 




0.001 0.01 0.1 

VJ 




0.001 0.01 0.1 



w 

Fig. 3.12. Vortex velocity / voltage power spectra and 
diffraction plots for increasing columnar defect lengths 
at a fixed vortex density p — 16/A — 0.00062/6q with 
(a) 1 = 1, (b) I = 3, and (c) I = 5. The results reveal 
how the system evolves from an ordered to a disordered 
configuration as the defect length increases. The wash- 
board signal decreases as the broadband noise grows. 



would repeatedly travel over the same pinning dis- 
tribution resulting in a periodicity corresponding 
to the 'time of flight' across the system. As a com- 
parison runs were performed for various point de- 
fect pinning strengths at the same low flux density 
p = 16 /A. The process was identical to the pre- 
vious section. However, due to the lower density, 
when the depth of the pinning potential well was 
increased to only 1.5 times its original value Uq, 
the washboard peak completely disappeared. The 
results of these simulations are shown in Fig. 13.131 
We observe that these results look quite similar to 
those from increasing the linear defect lengths. 

While the above power spectra are qualitatively 
comparable, the growth of the mean vortex line 
radius of gyration turns out to be quantitatively 
different in both situations, as shown in Fig. 13.141 
The data for columnar disorder initially grows as 
an exponential function, oc exp(— Zq/O where / is 
the length of the columnar defects, while the data 
for point defects is best described by a quadratic 
fit rg (X {U/Uq)'^. While we cannot offer a quanti- 
tative theory for these trends, they can be under- 
stood qualitatively by examining the average vortex 
velocity as these two pinning types are varied. For 
both cases as the 'pinning effectiveness' is increased 
the velocity decreases and the radius of gyration in- 
creases. We interpret the increase in radius of gy- 
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Fig. 3.13. Results for simulations with random point 
defects whose pinning strength is increased from its 
original value Uo in (a) to (b) U = 1.2 [/q and U — 
1.4(70 at a low vortex density p = 16/yl = 0.00062/6o. 



ration as the vortices being stretched as they are 
simultaneously held by the defects and pulled by 
the externally applied drive. As the linear defect 
length is increased the velocity decreases, but by 
length I — 20 (the height of the system) the vor- 
tices are still moving. [Compare the velocity trends 
in Fig. I3.14r a) and (b).] This indicates that the 
columnar defects cannot stop the vortices at this 
particular flux density and applied drive, and the 
'stretching' of the vortices is limited leading to satu- 
ration of the radius of gyration. On the other hand, 
for the point pins the pinning strength is not limited 
leading to greater stretching. 



4 Summary and Conclusion 

To summarize, a non-equilibrium Monte Carlo sim- 
ulation code has been developed to study the effects 
of different types of pinning centers on the velocity 
power spectrum of vortices driven through a disor- 
dered medium in the non-equilibrium steady state 
[23124] . Specifically, we have investigated the evolu- 
tion of the washboard signal as we have increased 
the vortex / magnetic flux density in the presence 
of point and columnar defects. In order to achieve a 
more complete understanding of the velocity power 
spectra, we have also examined the corresponding 
two-dimensional vortex structure factor and the av- 
erage flux line radius of gyration. 




1.1 1.2 1.3 1.4 1.5 1.6 

U/Uo 

Fig. 3.14. Mean flux line radius of gyration and aver- 
age velocity for (a) increasing linear defect length and 
(b) enhanced point pinning strength. The radius of gy- 
ration data for increasing lengths initially follows an 
exponential of the form rg oc exp(— Zo/0 while the data 
for increasing pinning strength is best described by a 
quadratic fit. A - average velocity, ■ - radius of gyra- 
tion {x component). 



We have confirmed that our numerical model 
displays the appropriate physical behavior over a 
large section of parameter space. For instance, vor- 
tices arrange into a six-fold lattice when interacting 
with only weak material defects. In the presence 
of an applied current and sufficiently strong pin- 
ning sites, the vortex lattice reorients with a lat- 
tice vector along the drive direction. Furthermore, 
we have obtained the current-voltage characteris- 
tics and found these to be qualitatively similar to 
TV curves obtained in experiments. In our simu- 
lation vortices depin above a critical applied force 
and gain velocity as the applied force is increased. 
Changes in the I-V curves due to different pinning 
site types and vortex densities occur as expected. 

At high applied drive, our algorithm produces 
saturating I-V curves, an artifact that originates 
in a necessary limitation of the maximum allowed 
displacements. However, we have provided evidence 
that suggests this limitation does not adversely af- 
fect the associated power spectra and other physical 
quantities in the regime investigated here. 
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Expected physical behavior is also observed in 
the vortex structure factor. For both pointlike and 
linear extended defects, as the vortex density is in- 
creased spatial ordering is observed to increase in 
the diffraction plots. When interacting with colum- 
nar pins, the vortex system structure factor is found 
to evolve from that of a liquid, to a smectic, and fi- 
nally a triangular lattice. For point defects only the 
regular triangular array is realized in the parame- 
ter space studied here. We remark that the struc- 
ture factor plots with point defects at low vortex 
densities appear qualitatively similar to those for 
columnar defects at higher flux densities. 

Velocity fluctuation or voltage noise power spec- 
tra measured parallel to the drive have been ob- 
tained for various vortex densities in the presence 
of both columnar and point defects. A narrowband 
signal is observed over a large vortex density inter- 
val, with the peak coinciding with the washboard 
frequency. The power spectrum has also been mea- 
sured perpendicular to the drive, and a signal at 
the washboard frequency is observed there as well. 
Harmonics have been detected at multiples of the 
washboard frequency for both types of disorder. 
For columnar defects the ratios of the power of the 
first and second harmonic to the third turn out to 
be larger than for point defects. However, as the 
density is increased, the ratio decreases for colum- 
nar defects. By varying the pinning strength of the 
point defects at a constant vortex density we have 
obtained similar harmonic ratios to that of colum- 
nar defects. This indicates that the harmonic ra- 
tios are not unique indicators of the type of mate- 
rial defects present in the sample. Rather, we think 
that the harmonic ratios should be a function of 
the degree of deformation of the vortex lattice. We 
remark that the detailed features of the flux flow 
power spectrum are also influenced by the config- 
uration of the leads above the sample surface 
In real experimental setups, these effects may mask 
the noise signatures observed in our simulations. 

In order to investigate the shape of the fluc- 
tuating vortex lines we have determined the aver- 
age radius of gyration in the presence of point and 
columnar defects in the x and y directions. The 
general (and expected) trend is for the radius of 
gyration to decrease as the density of the vortices 
is increased. Our results for the gyration radius also 
reveal that transverse fluctuations do not appear to 
play a large role in the thermal flux line wandering 
for the range of parameters investigated here: The 
transverse component in fact rarely exceeds the ra- 
dius of a pinning site except for the lowest vortex 
densities studied. 

In addition, we have measured these various ob- 
servables as the columnar defect length, i.e., the 
degree of correlation in the disorder was varied, 
and compared the results to the effects of varying 
just the pinning strength of point defects. Chang- 



ing the defect character from point-like to colum- 
nar shows similar results to increasing point defect 
pinning strength in both the diffraction and power 
spectral density plots. Different effects on the vor- 
tices through adjusting these two distinct pinning 
mechanisms become apparent upon comparing the 
growth of the radius of gyration in the direction 
along the drive. While the results seem to corre- 
spond to physical intuition, a precise theory ex- 
plaining the quantitative growth of the radius of 
gyration in either situation is currently lacking and 
will have to be developed. 
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